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Abstract 

Elaborating on a previous work by Marolf et al, we relate some exact results in quantum field theory 
. . . and statistical mechanics to the Bekenstein universal bound on entropy. Specifically, we consider the 

QP , relative entropy between the vacuum and another state, both reduced to a local region. We propose 

' that, with the adequate interpretation, the positivity of the relative entropy in this case constitutes a 

, well defined statement of the bound in flat space. We show that this version arises naturally from the 

original derivation of the bound from the generalized second law when quantum effects are taken into 
account. In this formulation the bound holds automatically, and in particular it does not suffer from 
the proliferation of the species problem. The results suggest that while the bound is relevant at the 
classical level, it does not introduce new physical constraints semiclassically. 
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1 Introduction 



The Bekenstein bound is a proposal for a universal bound on the entropy of a region in flat space. If S 
is the entropy, E is the total energy, and R is some typical size of the system, it writes 

S<\ER, (1) 

I where A is a numerical constant of order one. 
p<j ' The bound was originally conceived through a thought experiment involving black hole thermody- 

namics and classical physics [IJ. This is as follows. Consider an object much smaller than the black hole 
radius and drop it from a position outside the black hole. It will be swallowed and disappear behind the 
^ ^ horizon carrying away its entropy. The generalized second law of thermodynamics implies that the lost 
■ entropy have to be compensated (or exceeded) by an increase in the black hole area [2]. This is in turn 
OO . governed by the Einstein equations for gravity, and depends on the energy absorbed by the black hole, 
leading to ([I]). 

This derivation of the bound from the generalized second law of thermodynamics (GSL) was impugned 
^ . [U H] and defended [5l [6] in several occasions. However, despite this controversy around its association 
^ I with the GSL, the bound itself is interesting in its own right. The statement of the inequality ([T]) does not 
mention anything outside the area of ordinary flat space-time physics. After looking at it a bit closer one 
may convince oneself that what is behind protecting the bound are primarily the uncertainty relations 
(except for the proliferation of the species problem) [7]. This suggests that a full proof could be attempted 
in quantum field theory (QFT). 

However, the first and more difficult step in such direction is to properly define the quantities involved. 
In fact, the bound has been successfully tested in many different quantum systems, but the calculations 
have always involved at least some heuristic input. 

The main obstacle in properly stating the bound is the same that impeded the early attempts to 
construct a relativistic quantum theory for single particle (rather than a QFT). Eq. ([T]) requires defining 
a system which is bounded in space. This traduces quantum mechanically in the necessity of counting 
with a coordinate operator X. However, in a relativistic theory this would imply the existence of a time 
operator, a notion which clashes with quantum mechanics [8j. 
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Alternatively, one can try to define the localization of the system by imposing that the expectation 
values for local operators (i.e. the energy momentum tensor) in the region outside the one where the 
system is located coincide with the corresponding vacuum expectation values [9]. The entropy would 
in that case be given by its microcanonical expression S = log N(Eq) where N{Eq) is the dimension 
of the Hilbert space of states which fulfill these conditions, with energy below Eq. However, the set of 
states selected in this manner does not form a vector space, and the entropy cannot be computed |10j . 
Remarkably, the Reeh-Schlieder theorem implies that the vector space generated by this set of states by 
taking linear combinations span the whole set of states with E < Eq (including even very far located 
particles). The reason for this fact is that the states "localized" in this way are not really completely 
vanishing outside V. They have small tails outside V which can be exploited by adequately combining 
them (for a review on the problem of localization in QFT see [TT], for a different localization scheme see 
|12j . and for a review on the difficulties in correctly defining the bound see [13j). 

Therefore, we can compute the energy and the entropy for a global state, but in that case we do not 
know what to make of i? in ([T]). On the contrary, localizing the state in a QFT inevitably leads to pair 
creation, and the entropy and energy became ill defined. In a specific calculation this translates into 
divergences (think for example on the infinite entropy carried by the Unruh radiation in Rindler space). 
In order to eliminate them we can either impose a cutoff or ad hoc boundary conditions which eliminate 
the correlations of the localized state with the rest of the space. In both cases the result contains input 
from external elements and is ambiguous. 

A way out of this dilemma was discovered in [H] (see also p3]). There it was observed that a more 
reasonable measure for the entropy carried by the system in a region of the space is obtained by subtracting 
to the bare entropy of the local state the entropy corresponding to the vacuum fluctuations. This last 
one is present in absence of any energy, and is entirely due to the localization. The subtracted entropy 
turns out to be finite. 

A similar subtraction can be done to define a localized form of energy. Using this, we show how to 
define a version of the Bekenstein bound which has sense in quantum field theory. This coincides with the 
positivity of the relative entropy between two local states. We also show that this version of the bound 
emerges from known proofs of the GSL in the semiclassical quasistatic regime. 

The type of ideas involved in this work appears also in this context in the refs. [H [I6l [l71 [THl [19] , and 
especially in [T3], which largely motivated this paper. We hope it could help to clarify some aspects of 
the discussion in these works. 

2 A localized entropy 

First we focus on the left hand side of the bound inequality ([T]). Let us take a spatial region V lying on 
a Cauchy surface. In a quantum field theory there is an algebra of local operators ^(V) associated to 
V (or more precisely to its causal development). If p is the global state of the system, its localization 
pv to the region V is given by the restriction of p to A{V). In the continuum theory, it is not possible 
to give an entropy for pv since in general A{V) is a von Neumann algebra of a type which does not 
admit a trace [lOj . However, let us think in terms of a cutoff theory (i.e. with a lattice regularization), 
in such a way that the Hilbert space can be decomposed as a tensor product Ti. = Tiy ^ 'H-v of Hilbert 
spaces, associated to V and its complementary set —V on a Cauchy surface. These later give a product 
representation of the commuting algebras ^(1^) and A{—V). Then we have 

Pv = tr_y p , (2) 

where the trace is over Ti-v (that is, over the degrees of freedom lying outside V). In the limit of vanishing 
distance cutoff this state has a divergent entropy [20] 

S{pv) = -trpv- log Pv . (3) 

In order to see this, consider for example the global vacuum reduced to the Rindler wedge (where V is 
half a spatial hyperplane). S{pv) would be given naively by the entropy in a thermal gas at the Unruh 
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temperature T = {2it z)~^ , where z is the distance to the boundary plane. For massless particles the 
entropy density is then proportional to z~'^ and the entropy diverges proportionally to the integral of this 
quantity in half space, giving S ~ e~'^A, where A is the area on the boundary surface, and e is a short 
distance cutoff. 

This type of divergences is a general feature of the localized entropy. This is because the localization 
gives place to pair creation in the boundary region. For a general set V we have an expansion of form 

[2I1E21 

S{pv)=g2[dV]e-^ + gi[dV]e-^+g^[dV] \og{e) + {pv) , (4) 

where S^{pv) is a finite part, and the gt are local and extensive functions on the boundary dV, which are 
homogeneous of degree i. These divergent terms are ambiguous since they depend on the cutoff procedure 
(excepting the coefficient of the logarithmic term [21] 

To avoid the divergences one can put the system in a box to enforce boundary conditions and eliminate 
the entanglement with the exterior region. However, the box itself would be made of some material, and 
the consideration of the relevant fields for it would take us back to the same problem. Enforcing artificial 
boundary conditions cannot take away the ambiguities. 

However, the entropy one actually wants to measure is the one in the system under consideration on 
the region V and not the one due to the cloud of the vacuum fluctuations, which is present for any state. 
This is the crucial point made in [H]. Thus one should subtract from the system entropy the entropy in 
the vacuum state py = tr_y |0) (0|, corresponding to the same region V. It is important to note here 
that since the divergences have an ultraviolet origin the divergent terms in ^ are purely geometric and 
should not depend on the state p (at least if this later has a physically acceptable ultraviolet behavior, for 
example if it has finite energy). Thus, while both entropies contain divergences, the difference between 
the local entropies for two states is well defined, independent of the regularization. Our candidate for the 
left hand side of ([T]) is then 

Sv = S{pv) - S{p'y) . (5) 

When the bulk of the system is far from the boundaries of V the entropy Sy should approach the 
global state von Neumann entropy S{p). We think a formal proof of this could be attempted, for example 
considering the limit of sets with increasing sizes, or, for the Rindler wedge, the limit of large translations 
of the state p. A simple example was given in [Hj . where the entropy of a mixed state between different 
species of a single scalar particle in half space is computed. The result was there given in terms of a series. 
This is reviewed in the appendix, where we have also provided an analytic expression for the series sum. 



3 A localized energy 

Now we turn to the right hand side of ([T]). Again we find difficulties in trying to define the quantities 
involved. The first one is that, strictly speaking, the energy of a localized state like pv cannot be computed, 
since the Hamiltonian operator is defined in the whole Hilbert space, while the local density matrices act 
on the operator algebra generated by the observables localized in V. For this reason we need to define a 
local form of energy which has sense in the restricted space. 

Also, there is no general consensus on the correct value of A (27r from the original derivation [T], 
somewhat larger in [6]) or the meaning of R (circumscribing radius of the system in [1], the shorter 
dimension of the system in [23] ) . 

What we propose is in fact a candidate for the product \ ER in In this way we eliminate the 
ambiguities in the different quantities by a good definition of their produc10. 

In order to proceed we take the local vacuum density matrix (we still think in terms of a cutoff theory) 
and write it as 



e 



-K 



Pv = T—K ■ (6) 



tr e" 



similar step was carried out in [53] by using the harmonic resolution instead of the product E R. However, in the 
present case we do not want to abandon the localization of the states, which we consider an essential ingredient in the 
derivation of the bound from the GSL. 
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The hermitian and positive definite matrix py can always be written in this way since it has no zero 
eigenvalues (except when V is the whole space, where it becomes the vacuum state projector). The 
Hermitian and dimensionless operator K can be called the local Hamiltonian for V . This operator can 
also be defined directly in the continuum theory, where it is called the modular Hamiltonian [10]. Note 
that the operator e~^^'^ ~ (Py)*^ is unitary, and acts on the algebra of fields in 1/ as a kind of local 
dynamics with time s. 

When V is half a spatial plane (its causal development is the Rindler wedge) , the form of K is known 
exactly. In this case py can be written in the form ([6]) with K proportional to the boost operator which 
keeps the wedge invariant [25]. Explicitly, 

dxdy I dz z 7i{x,y, z,0) , (7) 

where x, y and z are the spatial coordinates, V is the set z > 0, and TC{x,y,z,t) is the Hamiltonian 
density operator. The local dynamics is then given by the trajectories of constant acceleration observers 

t' = cosh(27rs) t + sinh(27rs) z , 

z' = sinh(27rs) i + cosh(27rs) z . (8) 

Remarkably, this holds for any QFT. 

Our proposal for the product XER suggest itself from the formula ([7]) and is given by the expectation 
value of the operator K. When most of a system which has a size L is far from the boundary, R > L, 
this gives approximately K ~ 27rE R. 

As a different example for K consider the case of conformal field theories (in any dimensions), and 
when the set ^ is a sphere [26]. Then K is proportional to the generator of the conformal transformations 
which keep the sphere fixed. In polar coordinates it writes 

K = 271 dr [ r^dn^^^—^^^^^^n{r,e,(t>,0). (9) 







2R 

The local dynamics is shown in the figure 1 and is given by 



X' =R i(R + X±)-e'-^{R-X^)) 
± ^ ((fi + X±) + e2-(i?-X±)) ' ^'""^ 

where X+ = t + r and X_ = r — t are null coordinates, and here R means the radius of the sphere. 

Near the boundaries this flow approaches the one for the half space, and this feature is expected to 
be generic for any QFT and any set with smooth boundaries. Far from the boundaries the modular flow 
is not expected to be universal. For example, in the center of the sphere r = 0, at a distance R from the 
boundary, the velocity of the flow pO|) is half the one corresponding to a distance R from the boundary 
in the Rindler wedge. That is, one has K ~ ttER for a state well localized near the center of the sphere. 

For massive or interacting theories, or sets of more complex shape, the modular Hamiltonian is not 
expected to be local, and the local dynamics does not act as a simple space-time point transformation on 
the fields. For example, for a free massive Dirac fermion we have [27] 

K= f dxdy {x) log {Cy^ -l){x,y)^{y), (11) 
Jv 

where 

Cv(x,y) = /0|M/t(x)*(y)|0\ (12) 



is the field correlator inside the set V. This Hamiltonian is still quadratic, but it is in general non-local 
with respect to the fermion fields. 
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Figure 1: The horizontal interval represents a spatial sphere (only one spatial dimension shown), while 
the diamond shaped set is its domain of dependence. The vertical dashed curves are the trajectories of 
the point transformations induced by the local dynamics in a conformal field theory. The horizontal solid 
curves are the translation of the sphere for different internal times s (from bottom to top, the curves 
correspond to integer s from —5 to 5). 

A local Hamiltonian can also be defined for other states and space-times, in a way analogous to ([6]). 
One important case where it is possible to compute it exactly is when the state on V is thermal with 
respect to a given notion of time. In this case the reduced density matrix is of the form 



with E the time translation operator, and /3 the inverse temperature. Then, we have comparing with ([6]) 



This is the case for a black hole in the Hartle-Hawking state. For a Schwarzschild black hole we have in 
this case K = SttGM E, where E is the energy operator as measured by asymptotic observers, and M 
the black hole mass. More explicitly, in this case we have 



where Ti. = Tq is the Hamiltonian density in Schwarzschild coordinates, and the integration is over the 
exterior of the black hole region. As in the Rindler case, this is local and has a universal expression 
which is in principle valid for any theory. As it is well known, this operator approaches d?]), proportional 
to the boost operator in the Rindler space, very near the horizon. For this we have to pass to distance 
coordinates, dz = {1 — 2MG/r)~^^/'^^dr, and take the near horizon limit, (r — 2MG) ^ 2MG. In the 
opposite limit, r ^ 2MG, K is proportional to the usual Minkowski space time translations operator. 

Returning to flat space-time, and to the vacuum state induced local Hamiltonian there is one 
further problem with taking its expectation value on the localized states which we have not mentioned 
yet. To be concrete, consider the Rindler wedge. The vacuum state is thermal with respect to the boost 
operator with temperature (27r)~^. Thus, heuristically, for a massless field we can think in a thermal state 
for the usual time translations with temperature varying as T ~ l/z, and then "boost energy" density 
going as z^^. This diverges when integrated near the boundary. Another aspect of the same problem 
follows from the expression ([6]) for the density matrix in terms of K. The replacement K — > + constant 
does not change py. This makes K (and its expectation values) ambiguous by redefinitions by a constant 



Both of these problems can be cured by subtracting the corresponding vacuum expectation value of 
K in order to eliminate the divergent contribution coming from the vacuum fluctuations. The proposal 



K = (3E. 



(14) 




(15) 



term. 
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for the right hand side then reads 

Kv = ti{Kpv) - tr{Kpl) . (16) 

As happens for the entropy, this subtraction should lead again to a well defined cutoff independent 
quantity. 

4 Relative entropy and the Bekenstein bound 

Our version of the Bekenstein bound then reads Sy < Ky, which is 

S{py) - S{p^y) < tT{Kpy) - triKp^y) . (17) 

According to ([6]) we have K = — log py — log(tr e~^). Using this and tvpy = tvpy = 1, we can write pT|) 
simply as 

tv{py log py) - tr(py log p^) > . (18) 

This is just the statement of the positivity of the relative entropy S{py\py) between the local density 
matrices corresponding to the system's state and the vacuum, both reduced to V. 
The relative entropy between two density matrices pi and p2 is defined as 

S{pi\p2) = tr(pilogpi) -tr(pilogp2) , (19) 

and has the meaning of the statistical distance between the pair of stated [28j . It is always positive 

S{pi\p2)>0, (20) 

and is zero only when pi = P2- This gives a proof of the present formulation of the bound. As an example, 
for an object far enough from the boundaries in a half space the relation (|17p approaches ([T]) with A = 27r. 
In this sense (I17p captures what is underlying the validity of the bound in several examples studied in the 
literature. 

The relative entropy has also a well known application to physics in relation with the free energy. If 
the state p2 is thermal with respect to the Hamiltonian E as in (jl3p . its relative entropy with the state 
Pi is 

5(pi|p2)=/?(F(pi)-F(p2)), (21) 
where the free energy F(p) is given by 

Fip)=ti{pE)-TS{p). (22) 

Therefore, the positivity of the relative entropy is equivalent to the fact that the free energy of the thermal 
equilibrium state is the minimum one. 

One interesting point is how this formulation of the bound manages to avoid the species problem. 
The reason was first discovered in [T^j. The species problem in the Bekenstein bound arises since if the 
number of independent particle species M is increased we can make the entropy of py to increase keeping 
its energy fixed. One example is to take py = ^ py, where py represents a state containing only 
particles of the i type. The py for different i can be taken to be the same state with only the particles 
species interchanged. The entropy S{py) will then increase with M. However, the reduced vacuum 
state Py = ®py also becomes more entropic, and even if S{py) and S{py) both increase with M, their 
difference saturates for large enough values of the species number (a particular example of a single particle 

^In quantum field theory S{pv\pv) is perfectly well defined, independently of any cutoff procedure, even if the existence 
of the density matrices requires regularization. This mathematical definition is achieved by the Araki formula for the relative 
entropy for a pair of states in the local operator algebras [29]. According to the discussion in the previous sections the 
quantities Sv and Kv should also have a direct definition in the continuous QFT (perhaps it is even posible to define them 
for an arbitrary von Neumann algebra). However, to the author's knowledge these mathematical definitions have not been 
worked out in the literature. This interesting point deserves further study. 
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was studied in and is reviewed in the appendix). Moreover, the relative difference between the states, 
as measured by their relative entropy, diminishes with M, because it becomes more difficult to separate 
the state from the thermal atmosphere generated by the localization. To be more explicit, using the 
property of joint convexity 

S (E ^iPi I E ^^P2) < E ^^Sip\\p2) > (23) 
we have for the particular states considered above 

< S ^ pUp'v) < ^(pVIPy) • (24) 

The interpretation of these inequalities is that the relative entropy measures the degree of distinguishability 
between the two states [28], and, in contrast to the entropy, decreases under mixing, making the bound 
more strict. However, the bound never get violated since when the states are totally undistinguishable 
the relative entropy is zero. 

Another question is which region V makes the bound more severe for a given state. We have shown 
above that the sphere can make a better bound than half the space. In general we have that the relative 
entropy is monotonously increasing with size 

S (pvIp^v) < S {pwIp^w) ^ VQW. (25) 

Therefore its smallest values hold for small sets, where the state pv resembles the vacuum and S {pv\Pv) 
tends to zero. For large sets which include most of the state p it would increase without upper bound 
with the typical distance R between the set boundary and the system. Perhaps a more interesting and 
more difficult way to state this question is for which set the ratio Sy /Ky < 1 achieves a maximum for a 
given state. It is possible that in this case the set must adjust itself to some "form of the state", enclosing 
most of the entropy at the minimum distance. 



5 Discussion: the Bekenstein bound and the GSL 

The Bekenstein bound follows from the generalized second law by the thought experiment involving 
classical physics described in the introduction. However, there are, at the semiclassical level and in the 
quasistacionary regime, at least two different proofs of the GSL. These do not make use of a new principle 
like the Bekenstein bound, suggesting that properly including the quantum corrections the bound must 
hold necessarily. The purpose of this section is to make contact with these proofs. 

The Sorkin's proof for the semiclassical quasistatic case |17j relies only on the region of the space-time 
external to the black hole, which has autonomous evolution (this proof must not be confused with the 
more general proof-scheme of the GSL in quantum gravity by the same author given in [E]). Let pi 
and p2 be the two external states corresponding to the two different times t2 > ti, characterized by two 
Cauchy surfaces Si and S2. Sorkin also assumes the existence of the Hartle-Hawking thermal equilibrium 
state for the exterior region, /Ohh ~ e t ^ with T the black hole temperature and E the time translations 
symmetry operator for the static geometry. Here E/T is also the local Hamiltonian for the static exterior 
geometry in thermal equilibrium. Under the quasistatic approximation we can make the assumption that 
the time translation invariant state phh is the same on Si and S2. If the time evolution for the external 
states is completely positive-trace preserving (CPT) |28j we have 

'S'(piIphh) > S{p2\prr) ■ (26) 

In fact, the CPT property in this case holds, and coincides with the monotonicity for the operator algebras 
(f25|) . taking into account that the state p2 is obtained from pi by unitary evolution from ti to t2 plus 
trace over the degrees of freedom which on S2 lie inside the horizon. This inequality is equivalent to 

5(P2)-5(pi)-M2_Mi)>o, (27) 
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and using the quasistatic relation AS'bh = ((-£■) i — {^)2)l'^i gives the GSL. 

Now we turn to the experiment of letting an object be swallowed by the black hole. The least 
constraining inequality from (f26l) is just ^(piI/jhh) > 0, which happens for p2 = Phh- This gives our 
proposal (jl7p for the Bekenstein bound, but here the role of the vacuum is played by the Hartle-Hawking 
state. More explicitly, in this case ([27|) gives 



Both members of this inequality should be finite. As mentioned, the combination E/T is just the local 
Hamiltonian (jlSp for the external black hole geometry, and for objects near the horizon it approaches 
the corresponding one for the Rindler wedge. Thus, in the large black hole-flat space limit, where the 
space-time curvature can be neglected, we have the inequality (jl7p . where K is the local Hamiltonian of 
([7]) corresponding to the Rindler wedge. Therefore, the conclusion is that a form of the Bekenstein bound 
is indeed required by the GSL, but it is such that it does not introduce new physical constraints to QFT. 

The relation (j28p is a particular case of a more general one, where the Hartle-Hawking state is replaced 
by the thermal equilibrium state at temperature T ((see (I2ip ). The interpretation of the inequality in this 
case is just that any object may be made to come to thermal equilibrium with a thermal bath without 
violating the ordinary second law. 

However, in the thought experiment leading to the Bekenstein bound, the localization of the states 
introduces a very important difference: crossing the horizon means in this language to be absorbed by 
the thermal reservoir, and not just to be put into thermal contact with it. This may be interpreted as a 
sudden change in the Hamiltonian, rather than an interchange of energy, and in this case the result of the 
experiment cannot be explained thermodynamically. As the metric is the same for any type of matter, it 
may allow any object to be swallowed by the black hole (according to several works [31 H] this is a weak 
point of the Bekenstein deduction of the bound from the second law, since under certain conditions this 
could be impeded to occur). But at the same time, this same universal character of the metric gives place 
to a thermal atmosphere for any type of matter, which is already present alongside with the object initial 
state, and will also be present in the final state. The thermodynamic interpretation is then recovered by 
thinking in the autonomous exterior region alone, where an object and the thermal atmosphere together 
are seen to approach the state of thermal equilibrium. The Bekenstein bound then resolves in that for 
very entropic objects a very entropic thermal atmosphere is expected, and one has to take into account 
its entropy for the external region. The importance of thermal atmosphere then relies not on its possible 
mechanical effects on the object, but rather on the information available to the observer, which in some 
situations may not easily distinguish the object from the thermal cloud [13] . Note also that the thought 
experiment described above can be repeated word by word for any set V in flat space, where the exterior 
region is taken as the domain of dependence of V and /jhh is replaced by py. 

Another proof of the GSL has been given by Frolov and Page [18j . They consider the eternal black 
hole geometry for the quasistacionary process. Without entering into the details, this proof, apparently 
different in nature to the one by Sorkin, also boils down to the positivity of a particular relative entropy 



where here pi is the part of the incoming state which goes down to the future horizon, and /Ohh is the 
outgoing thermal state on the past horizon. The future and past horizon Hilbert spaces are mapped to 
each other by time reflection symmetry. Except for considering the horizon states instead of the exterior 
states, the conclusion is the same as above. 

A word of caution about the range of validity of these proofs (see also [6j). Strictly speaking, the 
relative entropy for two states in a given region can only be defined in a fixed geometry, and becomes 
ill defined once the backreaction is included. This is because two different states correspond in general 
to two different metrics, making the quantities in (jl9p uncertain (or UV divergent) [30]. Then, in both 
of these proofs the static limit is crucial, what means a fiat space limit for the Bekenstein bound. For a 
more general situation, but still in the semiclassical case, it is possible that we have to resort to express 
the quantities involved in terms of the mutual information [30] . 




(28) 



S{pi\p]m) > 



(29) 
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6 Appendix: a scalar particle in the Rindler wedge 



In this appendix we review a calculation of the subtracted entropy Sy by Marolf et al |14j . which is 
very instructive regarding the solution of the species problem for the Bekenstein bound. We also give an 
analytic expression for their final result for Sy, and express it in terms the relative entropy, demonstrating 
the validity of the approximations used in that work. 

Consider a state consisting in a superposition with equal weights of a single particle state corresponding 
to M different free scalar fields. We compute the relative entropy between this mixed state and the vacuum, 
both reduced to V, which is taken to be half a spatial hyperplane. 

To begin with, we recall that in this case the Hilbert space can be decomposed into a tensor product 
of the states in the left half space and on the right one, 7i = TIl ® Ti-R. Accordingly, the Minkowski 
vacuum state for each of the M fields decomposes as 



|0) =n^,fc(l-e-'^)^ ^e-^|iV,iV) . (30) 

N=0 

Here some discretization of the wave packets is assumed. The dimensionless variable uj is 2tt times the 
eigenvalue of the boost operator, while k represents the additional quantum numbers which define the 
wave packets. 

Following [14j we make the simplifying assumption that the single particle mode which we are con- 
sidering, in addition to have positive Minkowski energy, is well localized inside the right half space. This 
localization forces one to take to ^ 1. Then we can express the Minkowski annihilation operator for the 
mode to an exponentially small error as 

Vl - e-^ 

In this equation all the wave packet creation operators are normalized in standard way, [a, a^] = 1. 
This mode then has boost energy highly peaked around uj/{2tt). This expression neglects contributions 
proportional to aj^ and ai which are exponentially small in uj. We have retained the term proportional 
to in order to have a |0) =0. 

The problem then decomposes into the different frequencies w. A finite entropy for the reduced states 
require the presence of an infrared and an ultraviolet cutoff in order to make finite the sum over the 
modes. However, for the quantities we are interested in, which involve differences between the state and 
the vacuum, the only relevant contribution comes from the mode which is excited in the one particle state 
with respect to the vacuum. Then, for ease of notation, in what follows we omit all other modes from 
consideration. With this convention the single particle state is 

oo 

at|0) = (l-e-"') ^ ViVTle-"^ |iV,A^ + l) . (32) 

N=0 

Now consider M independent fields. The vacuum state is |0) = (8>i |0)j, and we take for the particle the 
mixed state M~^/^ a| |0). From here, tracing over the left Hilbert space, we obtain the reduced density 
matrices. Adopting a vector notation N = {Ni, Nm) and calling M = X^fcli we have 

oo 

p'v = il-e-n^'Y.e-'^^\Ni,...,NM){Nu...,NM\, (33) 



N=0 

1 ,„., 

PV ■ 



^ oo 

e-(l _ e-)M+i ^ ^g-^AT 1^^^ ^^^^ ^^^^ ^^^1 _ (34) 



M 

N=0 



Both of these are diagonal. The entropy of the vacuum corresponding to the mode in consideration is 

SiPv) = M (-^ - log (1 - e--)) , (35) 
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which is the one of a thermal ensemble of M independent oscillators. It also follows that 

S{pv) = - f; ^e-(l - e-)*^+We-^log ('^e-(l - e-)*^+iAAe-^') . (36) 

This is the expression found in [14J. It can be easily evaluated excepting for the term proportional 
to X]^=o-^^~'^"'^ 1°§ (■^)- order to obtain an analytic expression for this sum we use the integral 
representation for the logarithm 

7V>0 N>0 

Then, defining 

and deriving with respect to uj we have the differential equation 

dq(u), a) .... 1 , . 

- (a - a) + ^^^^ -1 = 0. (39) 

Taking into account the boundary condition lim^j^oo 

g{uj, a) 0, the solution is 

/•oo 

g{u;,a) = -e^"-^^^ e(i-")s'(l - (1 - e"^)"^) . (40) 



Using this result in (fSHj) and ([S7|) we obtain 
Sv = S(pv) - S(pI) ^ 

log(A/) - l„g(e" - 1) + ^ + du— (^(^^3^ j - Ij . (41) 

If we consider the number of species M fixed and take the limit of a particle far from boundary, such 
that Me'"" <C 1, we have 

Sv ~ log(M) + 0(Me~"') . (42) 

In this limit the finite local entropy Sy coincides with the global state entropy log(M). Prom (I4ip it follows 
that the entropy increases monotonically with M. If the behavior for large M were still logarithmical, 
at some point the Bekenstein bound would get violated. However, from ()41|) we see that the entropy 
saturates in the limit of large species number and fixed w, Me""" ^ 1, 

A similar calculation gives for the relative entropy between (jSSp and ()34p 

SiPvlp'v) = - ^^^Py) - ^^Pv)^ ■ (44) 



The positivity of the relative entropy then also shows that (|43p is an upper bound on the local entropy 
Sv no matter the number of species. This ensures the Bekenstein bound is not violated. 
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